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Abstract. 

An extensive simulation has been carried out to estimate the physical 
interpretation of dynamical factors such as < >, < E^Rh > in terms of high 

energy interaction features, concentrated in the present analysis on the average 
transverse momentum. 

It appears that the large enhancement observed for < EhRh > versus primary energy, 
suggesting in earliest analysis a significant rise of < pt > with energy, is only the 
result of the limited resolution of the detectors and remains in agreement with a wide 
range of models used in simulations. 
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1. Introduction 


Interest in the search for additional information in the hadron component of EAS at 
mountain altitude in terms of very high energy interactions Q has been spurred by the 
necessity to have a continuity of multiproduction features between present colliders at 
yi <1.8 TeV and the future LHC high energy physics at ^/s = 16 TeV {^/s represents 
here the centre of mass energy for the nucleon-nucleon collision). 

We review here the data obtained from calorimeters, or equivalent detectors, on 
high energy hadrons (longitudinal development and lateral spread) coupled with air 
shower arrays, where the primary energy is estimated from the electron and the muon 
components. 

From qualitative considerations, the most energetic hadrons {Eh > 1 TeV) can be 
supposed to conserve some significant properties correlated with the dynamics of the 
earliest interactions: for instance, the relation between the transverse momentum pt, 
the height h of the interaction, the hadron energy Eh and the radial distance Rh-, 


Pt = 


Rh Eh 


has been suggested as a possible measure of pt 


( 1 ) 


Following this line of inquiry, special attention has been given to the measurements 
of quantities like < EhRh >, < Eh >, < Rh > (where <> are the mean values per 
hadron - see |^, |^) and Uh, the number of hadrons above the energy threshold Eth- 
A convergent and significant increase of < EhRh > was observed in the different 
experiments i I § and the earliest interpretation was a tremendous rise of the 
transverse momentum < pt >, above -^/s = 5000 GeV, the minimum for pt (occuring 
from equation 1 when the hadrons come directly from the first interaction) being 
situated around 2 GeV/c 0. 

In more recent analysis measurements of hadron lateral distribution in Tien Shan 
calorimeter (which unfortunately have never been repeated |^) have been interpreted 
in terms of increasing role of high pt processes [^. The problem still remains open, 
since the number of reports on the subjest is small and “measurements of hadrons 
suffer more than other EAS components from systematic errors which have to be 
carefully estimated, using simulations of the detector response” [l^ . 

In parallel with our review of the experimental data, we have carried out EAS 
simulations with the program GORSIKA version 4.50 Bill, regarded as containing 
the major features of the present collider physics, in the energy range 10® - 10® GeV. 

Our analysis has concentrated on the hadron lateral spread to explain the 
specific effects on the different factors combined in < EhRh >, such as the sum 
of energy deposited, EEh, the hadron distance, Rh, and the hadron content nh (with 
consequences for the containment as well as the resolution of the detector). In order to 
appreciate the complexity of the hadron cascade mechanism, and how it can conserve 
or smear out the original interaction characteristics, the progression of our analysis is 
organized as follows: 

• physical dependence of interaction parameters (multiplicity, inelasticity, pt 
distribution....) 

• dependence on primary spectrum and composition 

• dependence on data acquisition and experimental resolution. 
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The survey of this dependence has focused on the largest detector employed, i.e. 
the Tien Shan calorimeter of 36 m^ Q. 

2. Simulation of the hadron lateral spread. 

2.1. Assumptions on pt and CORSIKA options. 

The angular emission of secondary hadrons depends on the ratio pt/pi ~ Pt/Eh, 
where pi is longitudinal momentum in laboratory system. Its evolution from CM 
frame to laboratory system is dominated, after Lorentz transformation, by the well 
known “collimation effect”, the emitted jet becoming narrower as the hadron energy 
increases. The multiplicity grows as In^ s, inelasticity changes slowly and pt rises 
regularly |Q, but in reasonable proportion (circle shape points in figure |^). 



Figured. Dependence of < pt > versus primary energy with the data interpretation 

of ref. [g from < E^Rh >• Iri the window options in CORSIKA adapted from ref. 

0 : 

- 1. HDPM + CERN UAl 85 (solid line) 

- 2. HDPM + CERN UAl Mimi 96 (dotted line) 

- 3. VENUS (dashed line). 

O accelerator and collider data compiled in [ |ll| 

□ < Ef^ Rh > interpretation from the Tien Shan experiment in . 

Other effects, such as the increase in pt for p-A collisions or A-A collisions (as 
described for the “Cronin effect” by Schmidt and Schukraft [^), or correlated with 
jet production, or the increase in pt for diffractive events [l^ are not very important 
or can be disregarded according to their respective cross sections. 

The correlation between < pt > and multiplicity (or central rapidity density 
dn/dy) cannot be neglected in the treatment of individual events more correlated 
with the semi-inclusive data; the correlation of the UAl experiment, described in 
iQ, is included in CORSIKA package (for the option HDPM) and has been used. 
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without change, in the main part of our simulation. We have, nevertheless, performed 
some simulations with a pt generator following the most recent correlation < pt > 
with dn/dy from the CERN UAl MIMI collaboration ||^, characterized by a greater 
rise of < Pt > at large multiplicity. In this way, we ascertained that the subsequent 
relative changes in < EhRh > are of second order when compared to the attempt to 
explain < EhRh > with < pt > values exceeding 2 GeV/c at ^/s = 5 TeV Q (squares 
on figure p. 



Figure 2. Dependence of < > versus Ne (for Et^ =1 TeV). 

- exact simulation (HDPM) output (dashed line) 

- after reconstruction described in section 3 (solid line) 

n Tien Shan data |l 7| 

O Tien Shan data 

Furthermore some simulations have been done with this extremal value 
< Pt > = 2 GeV/c around a primary energy Eq = 10^ GeV, after a special tuning 
of the Pt generator in GORSIKA: the options currently used have been HDPM in the 
usual form Q, as well as coupled to the pt generator elaborated following UAl MIMI 
results [^, and also the VENUS option |l^. A detailed description of the evolution 
with energy of the global interaction parameters used in the options of GORSIKA can 
be found in references - | p^ . 

The physical area explored following our selection of input parameters concerns 
primary energies between 10^ - 10® GeV, zenith angles between 0 and 30 deg for 
primary masses A = 1, 4, 14, 28, 56. The altitude for observation levels has 
been set to 3200 m with statistical samples of 500 simulated cascades per primary. 
Showers induced by heavy primaries have been simulated simultaneously with VENUS 
option and HDPM combined with superposition or , alternatively, with the abrasion- 
evaporation procedure [l^, describing the fragmentation of heavy projectiles. 
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2.2. Consequences for the hadron lateral spread. 

The results obtained for < EhRh >= {'EEhRh)/nh are shown in the figure || (dashed 
line), where < E^Rh > is plotted versus electron size Ng, for proton primaries and 
for a hadron energy threshold Eth = 1 TeV. The expected decrease versus size, as 
in previous calculations Q, is obtained. The agreement is tolerable only for the 
experimental data ||l^ at Ng near 2 • 10®. The lateral hadron densities A^, have been 
described by the experimenters with a simple exponential law: 


A/j = 


nh 

2ttRI 


Ho 


( 2 ) 


We have used the same functional dependence to fit the hadron distribution in each 
individual shower with two free parameters: normalization related to the total number 
of hadrons nt and mean radius Rq. The total number of hadrons < nu > varied from 
0.05 for HDPM model and primary Fe to 0.9 for VENUS option and primary proton 
for Ng = 10® and respectively from 3.5 to 9.0 for Ng = 10®. The average mean radius 
as a function of the shower size is presented in figure Similar remarks as in figure ^ 
are suggested by figure |^, where i?o decreases versus size, opposite to the general 
experimental tendency. 



Figure 3. Dependence of < Rq > versus Ne- (Lines and experimental data as in 
figure |). 
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Table 1. Dependence of < Ef^Rf^ > and < Rq > on interaction and shower 
parameters from the simulations as noted. 


Factors 

A 

< EhRh > 

< Ro > 

comment 



(TeV m) 

(m) 


a) primary mass A 

1 

3.16 

1.6 

HDPM -1- 

at Eo = 10® GeV 

4 

3.4 

1.70 

abrasion- 


14 

3.67 

2.21 

evaporation [|l6| 


28 

4.29 

2.52 

+ UAl MIMI 


56 

4.7 

3.24 

96 @ 


56 

3.81 

1.34 

HDPM+superp. 


1 

3.73 

1.09 

Venus 


Eth 

< EhRh > 

< Ro> 

comment 


(TeV) 

(TeV m) 

(m) 


b) hadron 

0.32 

1.988 

1.24 

HDPM 

threshold energy 

1.0 

3.34 

1.01 


(primary protons 

3.16 

5.51 

0.765 


of 10® GeV) 

10.0 

8.29 

0.57 



*0 

z 

< EhRh > 

< Ro> 


(gcm-^) 


(TeV m) 

(m) 

c) initial conditions, 

5 

2 

5.3 

3.08 

multiplicity 

5 

1/2 

5.6 

2.78 

(z = n^h/ < rich >) 

150 

2 

2.4 

0.98 

and depth (to) of first 

150 

1/2 

2.5 

1.26 

interaction (primary 





protons of 10® GeV) 


option: 

HDPM 



Eo 

< EhRh > 

< Ro> 

option 


(GeV) 

(TeV m) 

(m) 


d) primary energy 

10® 

4.248 

1.122 

HDPM 

(for protons 

310® 

3.678 

1.08 


and for E^fi = 1 TeV) 

10® 

3.58 

1.03 



310® 

3.34 

1.01 



10^ 

3.00 

0.95 



10® 

2.69 

0.93 



The total number of hadrons Uh reflects also a disagreement, characterized by an 
excess in the simulation [p0| , increasing with primary energy: the amplitude of this 
disagreement, when compared to the influence of primary energy and interaction pa¬ 
rameters, is illustrated by the values calculated in table 1. 

Setting < pt > above 2 GeV/c around a primary energy Eq = 10^ GeV gave no 
rise in < Eh Rh > and < Rq >: a slight decrease was even observed due to the larger 
transverse mass per secondary and to the consequent enhancement of the inelasticity 
in the former collision, reducing in proportion EEh at mountain level. 

Table 1 represents a small part of the total simulation which contains similar ef¬ 
fects to the examples presented. The total sample, up to 10^^ eV, confirms general 
aspects of the dependence, that we can summarize as follows: 
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• < i?o > and < EfiRh > decrease with energy even ii < pt > is rising with energy. 

• This effect persists even for extremal value of < pt > = 2 GeV/c and the 
“collimation effect” dominates over all the energy interval considered. 

• The height of the first interaction has an effect on the lateral spread, increasing 
< EhRh >. 

• The factors < Ro > and < EhRh > increase with the primary mass A. 

3. Event topology and detector response. 

The conditions of reception of the main hadron detectors employed at mountain alti¬ 
tude are compared in table 2. 


Table 2. Hadron detectors at mountain altitude . 


Location 

Detector 

Area 

(m2) 

Energy range 
(TeV) 

Resolution 
(cm ) 

Altitude 
(km a.s.l.) 

Comment 

Tien Shan 

ionization 

36 

> 0.3 

25 

3.2 

axis 


calorimeter 





contained 

Chacaltaya I 

em. chamber + 

8 


50 

5.2 



burst det. 






Chacaltaya II 

55 

8 

> 0.01 

50 

5.2 


Norikura 

proportional 

25 

0.02 - 0.12 

50 

2.77 

20m < r 


counters + r < 150m 

water tank + 
lead plates 


The hadron energy in the Tien Shan ionization calorimeter, composed of 
850 g cm“^ lead, is measured at a depth between 130 - 850 g cm“^ with an ac¬ 
curacy of 20% iQ. This thickness represents only 3.7 A/ (for an interaction length 
in lead A/ = 194 g cm“^) corresponding to a containment of about 80% of the cas¬ 
cades induced by individual hadrons [^. The coordinates of the single hadrons are 
measured with an accuracy of 25 cm in the 744 ionization chambers of the calorime¬ 
ter: those chambers, with an area of (25 x 300) cm^, originally arranged in 15 rows, 
are intended to measure the hadron number with an accuracy of 5% [^. Correction 
coefficients taking into account the limited lateral resolution and the possibility of ag¬ 
gregation of several hadrons within the chamber size bin (0.25 x 3) m were introduced 
by Danilova to recalculate the hadron energy spectra. The original method to 
separate and count the hadrons was based on the histograms of the ionization energy 
distribution in different rows in X and Y ||l^ . The limited efficiency of this treatment 
can be seen in figure which shows real contained target diagrams for hadrons above 
1 TeV and the same events reconstructed with the Tien Shan procedure; other diffi¬ 
culties appear also in the geometrical containment, the fraction of hadrons registered 
in one shower decreasing strongly with the primary energy (from 60% at Ne = 10® 
down to 20% at = 10® according to our simulation). 

The events simulated for section 2 have been treated by the numerical procedure 
of reconstruction of experimentally recorded detector response. The energy of each 
hadron was added to a respective ionization chamber output. As it is shown in 





DD 
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figure ^ from a detector plane four one-dimensional projections are obtained (each 
chamber is 3m long while the detector size is 6m x 6m). Then on each projection the 
possible “hadronic cluster” coordinates have been determined. Using as a separation 
criterion the dip in the projected energy deposit distribution of at least 70% of 
neighbouring chamber outputs, close chamber signals could be combined just as it has 
been performed in the Tien Shan experiment. Then after comparing the amplitudes 
on respective x- and y-projections taking into account the possibility of overlapping 
of some hadrons in one of the projections the so-called “hadronic clusters” have been 
found. They will be called hereafter the reconstructed hadrons. They are shown in 
the right plots in figure 

Our reconstruction procedure includes also hadrons of energies below the threshold 
(ITeV) which are mostly added to the large hadronic clusters. The correction 
coefficient of 2.64 (equal exactly to the one used by the Tien Shan group) have been 
applied to the cluster energy while converted to the reconstructed hadron energy. 
The position of the reconstructed hadron was the energy weighted mean of chamber 
outputs in each projection. As the shower axis coordinates we have used the most 
energetic reconstructed hadron coordinates just like it has been done in Tien Shan 
experiment. 

For all situations considered in section 2 the general decrease for < EhRh > and 

< > vs size turns to a large increase (solid line in figure ^ and figure as well as 

< EhRh > vs Eth (figure H). This is confirmed (figure ^) by comparison with the most 
recent data ||^ where the individual reconstruction indicates how the hadron content 
is underestimated near the axis. To preserve the clarity of the graphs, we have not 



10 1 10 10 " 


Eth [TeV] 


Figure 5. Dependence of < Eh Ru > on the energy threshold Eth- The squares are 
experimental points from M and [|L8[|. The energy dependence of the simulation is 
parallel to the experimental behaviour. 
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plotted on figures and ^ the results of the experiments §, 0, @-|| confirming 
the tendencies shown. Those values normalized to Tien Shan data can be found in 
reference 0 and reference 



Figure 6. Hadron lateral distribution function f(R) = dn(,/dR (for E}^ > ITeV) 
simulated and recalculated at 10® and 10^ GeV. Experimental data [Q. 


4. Discussion 

The proportion between the factor < Rq > calculated in table 1 and the resolution 
distance quoted in table 2 shows the role of underestimation of hadron densities near 
the axis (when several hadrons enter simultaneously in the same chamber bin). This 
effect will be even greater than in the Tien Shan ionization calorimeter in all the other 
detectors quoted . The values of Eth correspond to larger hadron densities near the 
axis, increasing the probability of aggregation in the same chamber bin. Also for larger 
resolution distance (i.e. wider bins) this probability increases. 

The simulation suggests that the hadron component and its lateral spread needs to 
be measured with more accurate calorimeters and a resolution distance substantially 
smaller than 25 cm. 

We have shown that the rise < EhRh > is a consequence only of the complex 
treatment of the data which is inadequate for the superposition of hadrons, and the 
rise of < pt > derived from < EhRh > and shown in the figure is not valid. The 
lateral spread of muons above 200 and 500 GeV, also correlated with < pt > in 
earliest interactions, doesn’t require such an enhancement of < pt > | |^ ; the main 
bulk of these lateral muon distributions being explained with a constant value of 
<pt> = 0.4 GeV/c II . 









11 


As follows from our simulations, even if < pt > is rising with energy, the hints 
in < iio > and < E^Rh > cannot be expected to reflect the pt behaviour, as they 
depend on several factors of the hadron cascade developement. 

To explain the enhancement in figure an extremely large increase in primary 
mass, proportional to the primary energy, was introduced in the simulation Q| to 
transform the general decrease of < EuRh > versus Ne in a visible increase. Such 
conjectures are no longer necessary and the agreement can be obtained without a large 
increase of the primary mass. 

The Pt dependence used here in the 3 options of CORSIKA covers a wide range of 
models; for instance, the values of pt in option 2 are very close to the values used in 
SYBILL as shown in ||^ by Knapp, Heck and Schatz, and we can conclude that 
the present measurements of hadrons at mountain altitudes are in agreement with 
a majority of models used in EAS simulations. They require neither a large < pt > 
enhancement nor a noticeable increase of primary mass. Their significance will remain 
limited until the resolution distance of the detectors is seriously improved. As it was 
advanced by Clay , “the results in hadron measurements depend critically on how 
well individual hadron signals are resolved in the hadron detector” and this has been 
confirmed by the present work. 
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